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ABSTRACT 


The  focus  of  the  study  during  the  last  three  years  was  to  analytically  understand  the 
effects  of  material  properties  as  well  as  configurations  on  the  magnetoelectric  coupling  in 
monolithic  mangetoelectric  materials  and  magnetostrictive/piezoelectric  layered 
composites  (denoted  as  MELC).  The  models  generated  during  the  last  three  years  provide 
valuable  information  on  understanding  the  magnetoelectric  coupling  behavior  and  hence 
maximized  the  magnetoelectric  coupling  for  magnetoelectric  system.  The  first  part  of  the 
efforts  had  been  devoted  to  developing  continuum  level  models  for  both  linear  and 
nonlinear  single  phase  magnetoelectric  materials,  including  the  solutions  to  representative 
problems.  For  the  second  effort,  we  focus  on  the  MELC  modeling.  A  total  of  six  MELC 
configurations  were  studied  by  this  model,  including  three  field  orientations,  longitudinal, 
transverse,  and  in-plane,  in  both  1  -D  and  2-D  plane  geometries.  By  using  the  modeling 
analysis,  a  3-D  a'  design  map  covering  the  span  of  compliance,  Poisson's  ratio,  and 
piezomagnetic  coefficient  ratio  of  the  magnetostrictive  phases  was  generated.  We 

also  manufactured  and  tested  several  MELC  systems  with  reasonable  agreement  betw'een 
experiment  and  theoretical  predictions. 


ACHIEVEMENTS 
Analytical  Modeling 
Monolithic  Magnetoelectric  Materials 

This  portion  of  the  research  focuses  on  developing  new  modeling  capabilities  for  a 
heterogeneous  layered  thin  film  magnetoelectric  and  coupling  these  solutions  with  a  new 
continuum  model  to  predict  the  dynamic  response.  The  first  model,  which  is  quasi-static, 
is  used  to  understand  the  complex  internal  interactions  intrinsic  to  monolithic 
magnetoelectric  materials. 

Initially  a  continuum  level  constitutive  model  was  formulated  for  multiferroic 
materials  using  an  energy  approach.  The  linear  model  was  used  in  the  governing  equations 
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and  a  proof  demonstrating  uniqueness  was  obtained,  The  general  formulation  was  applied 
to  several  standard  problems,  one  of  which  is  described  below. 

One  analysis  studied  a  rectangular  sample  of  multiferroic  material  occupying  a  region 
of  two-dimensional  space  denoted  by  0  <  jc  <  /  and  |  z  |  <  ^,  as  illustrated  in  Fig.  1  (b).  The 
material  has  been  polarized  and  magnetized  through  the  thickness,  i.e.,  parallel  to  the 
3-direction  or  z-a.xis,  and  is  assumed  to  be  transversely  isotropic  in  the  1-3  plane.  The  body 
is  assumed  to  be  in  plane  stress  with  constant  field  values  in  the  width  or  2-direction.  On  its 
boundaries,  the  material  occupying  the  two-dimensional  region  described  is  subjected  to 
uniform  electric  voltage  (±Fb)  as  well  as  a  uniform  magnetic  potentials  (±f^o)  at  its  top  and 
bottom  surfaces  (z  =  ±h),  respectively.  The  material  is  subjected  to  a  linearly  varying  stress 
((T,)  through  the  thickness  of  the  beam  at  the  end,  representing  an  applied  moment  at  its 
ends  (x  =  0  and  x  =  l). 


Figure  1.  (a)  Loading  Scenario  Considered  for  Analysis:  A  Two-Dimensional  Beam 
Experiencing  Combined  Loading,  (b)  Loading  Scenario  Diagram. 


Using  the  derived  constitutive  relations  and  the  boundary  value  problem  described 
above,  relations  describing  the  ability  of  the  material  to  convert  and  store  input  energy  or 
energies  for  the  3 1  -directional  case  was  computed.  This  came  in  the  form  of  an  effective 
coupling  coefficients  obtained  through  the  analysis,  given  in  Eq.  1  in  terms  of  both  material 
properties  (la)  and  the  corresponding  individual  coupling  coefficients  (lb)  for  bipartite 
interactions. 
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The  electro-magnetic  coupling  coefficient  is  of  particular  interest  in  this  study.  The 
square  of  the  effective  coupling  coefficient  (^31,8^)  when  plotted  as  a  function  of  the  square 
of  the  electro-magnetic  coupling  coefficient  (^3 !,£:«)  reveals  interesting  information  about 
the  coupling  behavior  of  the  material.  It  can  be  seen  in  Fig.  2  that  the  maximum  energy 
converted/stored  (i.e.,  maximum  )  by  the  material  system  occurs  at  values  of  /tsi.fw 

close  to  zero  and  one.  It  is  in  this  region  of  Fig.  2,  where  there  is  little  transfer  between 
electric  and  magnetic  energies,  and  as  a  result  the  coupling  exhibited  in  the  material  may  be 
approximated  by  a  superposition  of  piezoelectricity  and  piezomagnetism  results. 


^31, 
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Figure  2.  Effective  Coupling  Squared  Plotted  as  a  Function  of  Electro-magnetic  Coupling 
Squared. 


Magnetoclectric  laminate  composites  (MELC) 


In  the  program,  a  theoretical  model  was  also  developed  for  the  predicting  the 
magnetoelectric  (M-E)  voltage  coefficient  for  a  magnetostrictive/piezoelectric  layered 
composite  (MELC).  Three  field  orientations;  including  longitudinal,  transverse,  and 
in-plane,  with  two  in  plane  geometries,  2D  and  ID,  were  studied  with  respect  to  the  volume 
fraction  of  piezoelectric  phase.  As  shown  in  Fig. 3,  three  different  field  orientations  are 
studied  consisting  of  (i)  poling  direction  and  applied  magnetic  field  being  through  the 
thickness  (longitudinal  mode,  L);  (ii)  poling  direction  being  through  the  thickness  and 
magnetic  field  applied  in-plane  (transverse  mode.  T);  and  (iii)  poling  and  applying 
magnetic  field  in-plane  (in-plane  mode.  In).  For  each  of  these  modes,  ID  and  2D 
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geometries  are  studied;  thereby  representing  a  total  of  six  configurations.  In  discussing 
these  cases,  nomenclature  ofL2D,  LID,  T2D,  TID,  ln2D,  and  In  ID  are  used  to  denote  the 
orientation  (L,  T,  or  In)  and  the  geometry  (ID  or  2D).  When  a  magnetic  field  is  applied  to 
each  geometry,  a  deformation  is  generated  in  the  magnetostrictive  layer  which  is 
transmitted  to  the  piezoelectric  layer.  This  results  in  charge  accumulation  on  the  electrodes 
and  an  induced  electric  field  in  the  piezoelectric  layer.  The  magnetoelectric  voltage 
coefficient  a  defined  as  5E/5H  is  used  to  predict  magnetoelectric  coupling  effect  as 
presented  in  the  Fig.3. 
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Figure  3.  Representative  configurations  with  predictions  of  magnetoelectric  voltage 
coefficients  (a) 
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During  the  third  year,  a  wide  range  of  studies  were  conducted,  varying  the 
magnetostrictive  material  properties  and  the  volume  fraction  of  each  configuration.  This 
study  enabled  us  to  determine  the  individual  effects  of  MELC  configurations  and 
magnetostrictive  materia!  properties  on  the  homogenized  M-E  voltage  coefficient,  a', 
curves. 
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Table  I  Material  parameters  (compliance  s,  piezomagnetic  coefficient  q,  piezoelectric 
coefficient  d)  for  lead  zirconate  titanate  (PZT),  cobalt  ferrite  (CFO),  nickel  ferrite  (NFO), 
lanthanum  strontium  manganite  (LSMO),  and  composite  Terfenol-D  (C.TFD)  are  used  for 
calculation 


The  effective  properties  used  to  predict  the  M-E  coefficients  are  shown  in  Table  2. 
Figure  4  shows  two  plots,  one  for  a  CFO-PZT  composite  laminate  and  the  other  for  a 
NFO-PZT  laminate.  Both  figures  plot  the  homogenized  M-E  voltage  coefficient  a’  for  the 
six  cases  as  a  function  of  PZT  volume  fraction  '’v .  The  magnitude  of  a’  ranges  from  1 10.8 
mV/cm  Oe  to  1037.0  mV/cm  Oe  for  CFO-PZT  and  2A. 9  mV/cm  Oe  to  375.3  mV/cm  Oe  for 
NFO-PZT.  The  larger  a’  values  in  the  CFO-PZT  system  compared  to  NFO-PZT  are 
attributed  to  the  larger  piezomagnetic  coefficients.  While  the  CFO-PZT  system  has  larger 
a’,  the  peak  values  arise  at  the  same  '’v,  this  is  due  to  the  same  compliance  values  for  CFO 
and  NFO  these  two  magnetostrictive  materials. 


Figure  4.  Homogenized  M-E  voltage  coefficient  oc'  for  CFO-PZT  and  NFO-PZT  an  a 
function  of  volume  fraction  of piezoelectric  phase 


5 


Another  interesting  feature  at  these  graphs  is  the  ordering  of  a’.  That  is,  in  the 
CFO-PZT  system,  the  largest  to  smallest  a’  (i.e.  lnlD>ln2D>T2D>L2D>Tl  D>L1D) 
occurs  over  the  whole  '’v  range;  however,  in  NFO-PZT,  the  ordering  changes  with  ^v. 
For  '’v  less  than  0.53,  the  ordering  is  lnlD>In2D>T2D>L2D>TID>LID  while  for  '’v 
larger  than  0.53,  it  is  lnlD>In2D>T2D>T]  D>L2D>LI  D.  Since  the  two  magnetostrictive 
materials,  CFO  and  NFO,  have  the  same  compliance  values,  this  change  in  ordering  is 
attributed  to  piezomagnetic  coefficient. 

In  order  to  investigate  the  effect  of  piezomagnetic  coefficient  further,  both  the 
piezomagnetic  coefficient,  qj3,  but  the  ratio,  qjj/qji  is  considered.  Figure  5  shows  results 
for  a  CFO-PZT  system  with  two  different  qjs/qsi  ratios,  i.e.  -3.38  and  -2.0.  By  comparing 
these  two  plots,  we  can  see  the  largest  to  the  smallest  a*  for  each  case  changes.  For 
q33/q3i~3.38,  the  sequence  is  lnlD>ln2D>T2D>L2D>TlD>LiD,  while  for  q33/qji=-2.0 
the  sequence  is  ln2D>lnlD>L2D>T2D>TlD>LlD.  It  is  worth  noticing  the  largest  value 
a’  changes  from  ln2D  to  In  ID.  Furthermore,  L2D  dramatically  surpasses  T2D  for  the 
ratio  q33/q3i=-2.0  while  T2D  is  larger  than  L2D  for  the  ratio  q33/q3i=-3.38.  By  referencing 
to  Figure  1  and  equations  L2D  and  T2D,  we  conclude  that  increasing  the  ratio  q33/q3i 
significantly  improve  the  M-E  coupling  for  L2D  case  but  counteracts  T2D  response.  The 
result  is  also  in  agreement  with  the  M-E  effect  measurement  in  L2D  type  sample  studied  by 
Ryu  which  reaches  the  highest  M-E  coupling  (4.68x10^  mV/cm  Oe)  for  the  L2D  when 
compared  to  the  T2D  configurations. 


Figure  5.  Comparison  of  homogenized  M-E  voltage  coefficient  a'  for  CFO-PZT  under  the 
different  ratio  qifqai-S.SS  v.s.-2.0) 


On  the  contrary,  if  the  q33/q3i  ratio  is  fixed  while  only  q33  value  is  changed,  we  can 
study  the  effects  of  q33  on  a’.  As  shown  in  Fig.  6,  by  increasing  the  q33  of  NFO  from 
-680x  [0"'^  m/A  to  -800><  1 0'*^  m/A,  the  peaks  of  each  configurations  increase  roughly  by  a 
factor  of  1.15  without  changing  sequence.  Therefore,  given  the  same  compliance 
properties,  the  ordering  of  a’  for  each  case  is  influenced  by  the  ratio  q33/q3i  instead  of  the 
q33  value. 
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Figure  6.  Comparison  of  homogenized  M-E  voltage  coefficient  a' for  NFO-PZT  under  the 
different  piezomognetic  coefficient  qa  (-680'>^I0'‘"  m/A  v.s.  -800^ I m/A) 

To  study  the  influence  of  compliance  property,  comparison  of  LSMO  and  composite 
Terfenol-D  (Terfenol-D  particulate  mixed  with  resin)  are  presented  in  Fig.  7.  LSMO-PZT 
exhibits  a  low  a’  ranging  from  1 6.2  mV/cm  Oe  to  90. 1  mV/cm  Oe  while  the  Terfenol-D 
composite  a’  ranges  from  73.0  mV/cm  Oe  to  676.4  mV/cm  Oe, 


700 
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Vofum*  fractton 
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Figure.  1  Homogenized  M-E  voltage  coefficient  a’  for  LSMO-PZT  and  Composite 
Terfenol-D-PZT  (C.TFD-PZT)  as  a  function  of  volume  fraction  of  piezoelectric  phase 


By  comparing  the  peak  values  of  L2D  in  LSMO-PZT  to  NFO-PZT,  we  can  see 
LSMO-PZT  is  lower  by  a  factor  of  2  (6 1.24  mP7cm  Oe  for  LSMO-PZT  and  102  mV/cm  Oe 
for  NFO-PZT).  However,  as  shown  in  Tablet,  the  magnetostrictive  property  qji  are  the 
same  for  LSMO  and  NFO  but  different  compliance  (LSMO’s  is  about  2  times  higher  than 
NFO’s).  Because  ot’  of  L2D  is  affected  solely  by  q^i  as  shown  in  Fig. 3,  the  weak  a’  is 
attributed  to  high  compliance  and  not  piezomagnetic  coefficient.  Further  comparison  is 
obtained  from  the  composite  Terfenol-D-PZT.  With  about  25  times  higher  piezomagnetic 
coefficient  (gjj=6,389xl0‘‘^  m/A  for  composite  Terfenol-D)  compared  to  qsi  of  LSMO 
(250><I0’'^  m/A),  the  a’  peak  of  InlD  case  (solely  affected  by  qjj)  only  exhibits  8  times 
larger  magnitude.  We  conclude  that,  as  the  compliance  further  increases,  not  only  the 
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magnitude  of  a’  decreases  but  a’  peaks  shift  to  smaller  values  of  ^v.  This  suggests  that 
mechanical  impedance  matching  is  an  important  property  to  maximize  a’. 

Results  show  that  the  piezoelectric  phase  volume  fraction  to  achieve  maximum  M-E 
coupling  effect  is  dependent  upon  the  compliance  of  magnetostrictive  phase.  Higher 
compliance  values  shift  a’  peaks  to  lower  piezoelectric  volume  fraction  and  attenuate  a’ 
values.  This  study  also  investigated  the  influence  of  the  piezomagnet ic  coefficient  qjj  and 
the  ratio  qii/qji  on  a’.  For  constant  ratio  of  qjs/qsi,  larger  qss  values  increase  the  M-E 
coupling  effect.  However,  changing  this  ratio  alters  the  relative  ordering  of  the  a’  values 
for  each  of  the  six  MELC  configurations  studied.  This  demonstrates  that  the  ratio  q^j/qji 
strongly  influences  the  selection  of  MELC  configurations  to  produce  the  largest  M-E 
coupling  effect. 

Using  the  results  from  this  modeling  analysis,  a  3-D  a'  sequencing  map,  as  shown  in 
Fig.  8,  covering  the  span  of  compliance,  Poisson’s  ratio,  and  piezomagnetic  coefficient 
ratio  qii/qit  of  the  magnetostrictive  phase  can  be  generated  to  aid  in  the  designing  MELC 
system.  Each  surface  presented  in  the  map  represents  a  transition  boundary  from  one 
sequence  to  another,  i.e.  a  transition  plane  on  which  a’  values  of  two  configurations  are 
equivalent.  For  the  region  between  surfaces,  we  have  provided  the  sequencing  that  is 
produced  for  each  particular  region.  The  order  of  the  a'  values  from  the  largest  to  the 
smallest  is  hence  clearly  illustrated  in  the  magnetostrictive  material  property  space.  As  can 
be  seen,  the  MELC  configuration  with  the  maximum  at'  value  changes  as  the  properties  of 
the  material  changes.  Therefore,  it  is  incorrect  to  assume  a  given  configuration  always 
exhibits  the  highest  a'  value.  Fig.  8  can  also  be  used  by  researchers,  manufacturers,  and 
designers  to  select  geometries  and  field  orientation  for  MELC  systems  that  achieve  the 
highest  magnetoelectric  effect  based  on  the  choice  of  magnetostrictive  phase  properties. 
Furthermore,  if  we  consider  the  variety  of  magnetostrictive  material  properties  in  MEMS 
fabrication  process,  this  mapping  can  also  be  employed  to  determine  expected  response 
from  a  wide  range  of  material  selection  and  process  setting. 
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Figure.  8  S-D  a'  sequencing  map  covering  the  span  of  compliance  '"^n ,  Poisson 's  ratio 
v^,  and  piezomagnetic  coefficient  ratio  of  the  magnetostrictive  phases 


Fabrication  and  Experiment 

A  Terfenol-D  layer  used  in  the  magneto-electric  laminate  composite  consisted  of 
aligned  Terfenol-D  particulates  and  a  polymer  resin.  The  fabrication  process  for  the 
Terfenol-D  layer  is  shown  in  Fig.  9.  The  1-3  Terfenol-D  composites  were  fabricated  by 
mixing  55%  by  volume  ball-milled,  50-300  pm  in  size,  Terfenol-D  particulates 
(ETREMA  Products,  Inc.,  Ames.  lA)  with  a  Vinyl  ester  resin  (Dow  Derakane  41 1-C50)  in 
a  Plexiglas  mode  with  a  cavity  of  24.45x1 4.45x12.55  mm^  The  mixtures  were  degassed  to 
eliminate  the  voids  inside  the  composite  and  the  liquid  slurry  was  placed  between  a  pair  of 
Nd-Fe-B  permanent  magnets  producing  a  magnetic  field  of  150  kA/m  to  align  the 
particulates.  The  composite  was  then  placed  in  a  vacuum  oven  at  70  “C  for  8  hours  to  fully 
cure  the  Vinyl  ester  resin.  Fig.  9(a)  shows  the  Terfenol-D  composite  after  demolding. 
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Figure  9.  (Color  online)  1-3  Terfenol-D  and  Magneto-Electric  samples  fabrication 


The  Terfenol-D  composite  shown  in  Fig.  9(b)  was  characterized  to  obtain  material 
properties  including  elastic  moduli  (“Tj/,  '"T;;),  piezomagnetic  coefficients  C^ih  "‘^12), 
and  relative  permeability  for  various  mechanical  loadings  as  a  function  of  a  bias 

magnetic  field.  The  data  is  important  to  ensure  that  variations  observed  between  test  data 
and  theory  are  not  the  result  of  inaccurate  properties  used  in  the  models.  Two  different  sets 
of  measurements  were  conducted  at  room  temperature.  The  first  set  of  tests  involved 
applying  a  bias  magnetic  field  (Hbias)  and  a  compressive  bias  mechanical  preload  (cri„oj)  in  a 
presence  of  a  time- varying  mechanical  loading  (u/ac)-  A  sinusoidal  <Ti,ac  of  2  MPa 
amplitude  was  applied  at  0.3  Hz  to  characterize  the  elastic  modulus  for  atias  combinations 
ranging  from  2  MPa  to  12  MPa  in  2  MPa  increments  with  Hbias  ranging  from  0  Oe  to  4000 
Oe  in  500  Oe  increments.  The  second  set  of  tests  applies  a  compressive  vvith  a 
time-varying  Ht,ae.  A  <7hms  ranging  from  0  MPa  to  12  MPa  in  2  MPa  increments  was 
applied  followed  by  a  sinusoidal  H/^gc  with  amplitude  of  4000  Oe  at  0.3  Hz.  For  both  sets 
of  tests,  data  was  recorded  for  7  cycles  and  average  material  properties  were  calculated 
with  trivial  statistical  deviation. 

Following  characterization,  the  composite  was  sliced  into  laminates  of  three 
thicknesses  of  0.69  mm.  1.34mm,  and  2.62  mm  with  the  same  24.45x12.55  mm^  planar 
area  as  shown  in  Fig.  9(c).  The  piezoelectric  ceramic  in  the  magnetoelectric  (M-E) 
laminate  composite  is  a  PZT-5H  plate  of  0.27  mmx24.5  mmxl2.6  mm  (Piezo  Systems, 
Inc.,  Cambridge,  MA).  As  illustrated  in  Fig,  9(d),  each  M-E  laminate  was  fabricated  by 
sandwiching  a  single  layer  of  Terfenol-D  composite  laminate  between  two  oppositely 
poled  PZT-5H  plates  with  a  76  pm  thick  silver  conductive  tape  for  bonding.  These  three 
M-E  laminates  were  cured  at  125  “C  for  2  hours  under  a  30  psi  compression  resulting  in 
M-E  laminates  of  various  piezoelectric  volume  fractions,  '’v ,  of  0.44, 0.29,  and  0. 1 7.  The 
piezoelectric  layers  were  electrically  connected  in  parallel,  Fig.  9  (e),  to  monitor  the 
voltage/charge  produced. 

We  fabricated  a  number  of  symmetric  MELCs.  Fig.  10  shows  a  sample  of  the 
fabricated  laminates  for  T2D  in  Figure  1 .  The  magnetoelectric  laminates  were  fabricated 
by  sandwiching  a  layer  of  1-3  Terfenol-D  particulate  composite  in  between  two  layers  of 
PZT-5H  to  eliminate  the  bending  of  the  MELC.  The  layers  were  bonded  using  a 
conductive  tape  with  a  thickness  of  25  -  50  pm  (Emerson  &  Cuming  CF3350).  The  MELC 
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were  cured  at  137  “C  for  1  hour  under  a  50  psi  compression  to  ensure  good  mechanical 
bonding  between  layers.  The  Terfenol-D  layers  of  different  thickness,  a  total  five 
piezoelectric  volume  fractions,  defined  as  the  thickness  ratio  of  PZT-5H  to  the  MELC, 
were  obtained  with  ’’v  of  0.17, 0.28, 0.34, 0.38,  and  0.44.  To  conduct  high  frequency  tests, 
the  magnetostrictive  layers  were  fabricated  with  Terfenol-D  particulates  of  50  pm  ~300 
pm  in  size  embedded  into  a  vinyl  ester  resin  matrix  (Dow  Derekane  41 1-C50). 


Figure  10.  Illustration  and  pictures  depicting  the  different  magnetoelectric  laminates 
manufactured 


Fig,  1 1  compares  the  theoretical  and  experimental  magnetoelectric  voltage  coefficient 
(a)  for  MELC  systems  with  the  volume  fraction  of  0.2,  0.4,  and  0.8.  The  solid  symbols 
mark  the  experimental  results  and  the  hollow  symbols  with  dashed  lines  represent  the 
theoretical  predictions.  All  the  a  curves  for  each  sample  show  similar  trend  which 
increase  to  the  peak  value  and  decrease  with  the  bias  magnetic  field.  Furthermore,  for  the 
case  of  ^v=0.2,  the  theoretical  and  experimental  a's  reach  the  identical  values  of  around  1 
V/cm-Oe  with  the  error  less  than  2  %.  This  provides  researchers  a  viable  method  to 
estimate  the  experimental  magnetoelectric  responses  by  using  theoretical  prediction. 
Although  the  results  are  promising,  there  is  a  deviation  for  the  other  cases.  As  the 
piezoelectric  volume  fraction  increases,  the  differences  between  the  test  and  theory 
become  more  significant.  This  complex  but  interesting  phenomenon  is  proposed  to  be 
caused  by  the  stress  transmission  lag  due  to  the  coupling  magnetostricti ve/piezoelectric 
interface  and  the  demagnetization  of  the  magnetostrictive  layer. 
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Figure  11.  The  theoretical  and  experimental  magneotelectric  voltage  coefficient  (a)  for 
magnetoelectric  laminate  composites  (MELC)  with  piezoelectric  volume  fraction  of  0.2. 
0.4,  and  0.8.  Solid  and  hollow  symbols  represent  test  and  theory  results,  respectively. 
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